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, We examine spherical gravitational collapse of a matter model with vanishing radial pres- 

' sure and non-zero tangential pressure. It is seen analytically that the collapsing cloud either 

(N 

, ^ ' tial density, tangential pressure and velocity profile of the cloud. A threshold of black hole 



formation is observed near which a scaling relation is obtained for the mass of black hole, 
assuming initial profiles to be smooth. The similarities in the behaviour of this model at 
the onset of black hole formation with that of numerical critical behaviour in other collapse 



(N 

^ ' models are indicated 
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O : I. INTRODUCTION 

u 

Astrophysical black holes have a minimum mass, known as the Chandrashekhar mass. Reahstic 
- ^ , matter in a star is almost stationary at the initial epoch of its collapse and has a characteristic scale 

^ ' which depends on the properties of matter. If one does not restrict to stationary realistic matter, 

then in principle one should be able to produce arbitrary small mass black holes from initially 
collapsing cloud, as there will not be any characteristic scale present in the cloud and the positive 
pressure in the cloud has tendency to push the matter in radially outward direction in spherically 
symmetric situation. 

There can be very many possible equation of states with pressure present that the collapsing 
matter can satisfy, but the non-linear Einstein equations are tractable analytically only for a very 
few equations of state. Here, we examine a spherically symmetric collapse model, in which the 
radial pressure vanishes and non-zero tangential pressure is present in the cloud. The Einstein 
cluster is such an example with vanishing radial pressure, where a cloud of rotating particles whose 
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motion is sustained by an angular momentum has an average effect of creating a non-zero tangential 
stress. Such a static system was first introduced by Einstein which was later generalized to 
non-static case [2]. Subsequently, the gravitational collapse models with a vanishing radial pressure 
have been investigated extensively by Magli 3| and others [J] for the purpose of examining the 
validity or otherwise of the cosmic censorship conjectures and to understand the black hole and 
naked singularity formation in gravitational collapse. 

We observe an intriguing feature in this model that it has a threshold of black hole formation 
and it shows a power law behaviour for the mass of the black hole near this threshold. We work in 
the comoving coordinates (t, r, 9, (p) and construct an effective potential for the collapsing shells, 
with the aid of which one can see transparently the various possible evolutions of the collapsing 
shells. We demonstrate for a set of initial data that the collapsing cloud either forms a black hole 
or completely disperses, depending on the values of initial parameters. We point out similarities 
in behavior of this model at the threshold of black hole formation with already established critical 
behavior seen in other matter models ^]. 

The outline of the paper is as follows. In Section II, we discuss the collapse equations and 
regularity conditions and in Section III, a tangential pressure model is constructed. In Section IV, 
various possible dynamical evolutions are illustrated and a scaling relation for the black hole mass 
is obtained. Discussion and conclusions are outlined in Section V. 

II. EINSTEIN EQUATIONS, REGULARITY AND ENERGY CONDITIONS 

We consider four dimensional spherically symmetric metric in comoving coordinates, 

where dCl'^ is the line element on two-sphere. The energy-momentum tensor is then diagonal for 
any general Type I collapsing matter field and is given as T/ = —p; = pr] Tg = = pa. 
This is a general class of matter fields that includes many known physical forms of matter [6|]. 
The quantities p, pr and pq are the density, radial and tangential pressures respectively. We take 
the matter field to satisfy the weak energy condition^ that is, the energy density measured by 
any local observer be non-negative. Then for any timelike vector we have Ti^V^V^ > 0, i.e. 
p > 0; p+Pr > 0; p+pe > 0. The dynamical evolution of the system is determined by the Einstein 
equations and for metric ([I]) these are given as 

F' F 
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R' 

v'{p+Pr) = 2{p0 -Pr)^ -P'r, 



2R> + R'^ + R^ = 0, 



G-H =1 



F 
R' 



(3) 
(4) 
(5) 



where (') and (') are partial derivatives with respect to t and r respectively and 

G{t,r) = e-^^R'^ H{t,r) = e-^" R^ . (6) 
The function F{t, r) is twice the Misner-Sharp mass for the collapsing cloud, which gives the 



a. 



Regularity at the initial epoch implies 



total mass within the shell of comoving radius r at time t 
F{ti,0) = 0, that is, the mass function vanishes at the center of the cloud. It is seen from ([2]) that 
the density of matter blows up when R = or R' = 0. Here R' = corresponds to a shell-crossing 
singularity. We use the scaling independence of the coordinate r to write, R{t, r) = rv{t, r) and we 
have 

v{U,r) = 1; v{tsir),r) = 0; v{ti,r) < 0, (7) 

where ti and tg are the initial and singular epochs respectively. The condition v{ti,r) < signifies 
initially collapsing shells. At the initial epoch we have R = r, i.e. v{t, r) = 1, and at the singularity 
R = and v = 0. At all other epochs v has a non-zero finite value for all r with < r < rb, where 
Tb is the boundary of the cloud. 

From the point of view of dynamic evolution of initial data from t = ti,we now have five arbitrary 
functions, given by i^{ti,r) = i^o{r),tp{ti,r) = ipoir), p{ti,r) = po{r),pr{ti,r) = Pro{r),pe{ti,r) = 
P0^^{r). They are all not independent, ([3|) gives a relation for them. To preserve the regularity 
and smoothness of initial data we assume that the gradients of pressures vanish at the center, 
P'roi^) ~ P'doi^) ~ ^- Thus we have a total of five field equations with seven unknowns p, pr, pe, ip, 
V, R, and F, giving a freedom of choice of two free functions to complete the system. This choice, 
subject to the given initial data and weak energy condition, determines the matter distribution and 
metric of the space-time, leading to a particular dynamical collapse evolution of the initial data. 
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III. TANGENTIAL PRESSURE MODEL 

In this framework, we now study a class of collapse models with vanishing radial pressure and 
non- vanishing tangential pressure, which have been studied extensively in recent years 0,0]. For 
the sake of transparency, we construct and consider below an explicit collapse solution with a 
non- vanishing tangential pressure. We choose the allowed two free functions, F{t,r) and u{t,r), 
in the following way. The vanishing or a constant pr implies that the mass function F has to be 
necessarily of the form F{t,r) = r^A^(r), where Ai can be any general function of r, and we also 
choose ^{t^r) = v^lR). 

Consider M now be of the form 

M(r) = moe-™2"', (8) 

where mo and m2 are positive constants. Using the above form of mass function in equation ([2]), 
we get Pr = 0, i.e. the radial pressure vanishes identically, and the density at initial epoch is given 
by po(ti,r) = (3 — 2m2r^)moe~'"^^^. We also take the initial tangential pressure to be of the form 

P0{ti,r)=ePO2'-' -1, (9) 

such that both pg and p'g vanish at the center. We note that while the above profiles are chosen 
in order to give one particular solution, the following analysis is valid for any smooth initial data 
with initial density and pressure expressed in even powers of r. In general, as f ^ 0, p ^ oo and 
the density blows up at = 0, which is a curvature singularity as expected. Using u = uq^R) in 
equation we get 

G(t,r) = /(r)e2-»(^), (10) 

where /(r) is another arbitrary function of r. In analogy with dust collapse models, we write 

f{r) = l + r^b{r), (11) 

where b(r) is the energy distribution function for the collapsing shells. We take it to be a smooth 
function, b{r) = — 6o + 62r^ + - • • . Using i/ in ([3]), the equation of state turns out to be 2pg = Rvfi p. 
Finally, using (flO]) in equation ([5]), we get 

/Ri? = _e''o(-f?V/(r)i?e2'^o -R + r^M. (12) 
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IV. DYNAMICAL EVOLUTION OF THE COLLAPSING SHELLS 

The initial density profile, tangential pressure profile and velocity profile of the cloud are now 
as specified above and we have to evolve the initial data to investigate the possible outcomes of 
the collapse. Initially, all the shells have the scale factor v{ti,r) as unity, with v{ti,r) < 0, i.e. an 
initially collapsing cloud. The possible bounce of a shell is given by the change in sign of v. 
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FIG. 1: Typical effective potential for a bouncing shell. 



The evolution of a particular shell is deduced from equation ([12]) • Rewriting ([12]) in terms of v, 
we get 

vi,^ = e^^^o [v j{r,v)+M] = -V{r,v), (13) 

where j{r,v) = [/(r)e2'^o(™) - Ij/r^. We caU V{r, v) as the effective potential for a shell. A similar 
effective potential has been constructed for studying critical behaviour in scalar field collapse in 
higher dimensions by Frolov [8]. As seen in FIG.l, the allowed regions of motion correspond to 
V{r, v) < 0, and dynamics of a shell can be studied by finding the turning points. Starting from 
an initially collapsing state {v < 0), there is a rebounce if -u = (at u = f 2 in FIG.l) before the 
shell has become singular. This can happen when V{r, v) = 0. Hence, to study the various possible 
evolutions for a particular shell, we need to analyze roots of the equation V{r, v) = 0, keeping 
the value of r fixed. We have confirmed numerically that there is no shell-crossing in the cloud 
for a wide range of initial data of interest, which suggests that if a particular shell with comoving 
coordinate bounces then all the other shells with r > must also bounce. This implies that if 
central shell at r = bounces, then the whole cloud must also bounce off. This suggests that to 
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investigate the situation when the whole cloud is just about to disperse off, it is sufficient to study 
the dispersal of the shells near the center. Therefore, to find the threshold of black hole formation 
and to get the scaling relation for the mass of black hole, we need to analyze the model only near 
the center. 

With the form of smooth initial data ([8]), ([9]), we can integrate ([3]) on t = ti, close to the center, 
and i^Q = uoiR) gives 



We neglected higher order terms above since we want to consider the evolution of shells near r = 
only. Near the center, equation (|13p can be written as 



where A{r,v) = 1 + (2p6i2r^f^)/3mo and P5 = {4:p0^ + mopg^lSpg^ + 10m2))/18mg. The first factor 
in V{r,v) i.e function A(r,v), being the \goo\ term, is always positive and does not contribute to 
the bounce of shells. The main features of evolution of cloud derive from the second factor in 
V{r,v), a quintic polynomial having five roots in general. Only positive real roots correspond to 
physical cases. Since V{r, 0) = —M. < 0, any region between R = rv = and the first positive 
zero of V{r, v) always becomes singular during collapse. The region between the unique positive 
roots is forbidden as < there. For a particular shell to bounce it must therefore lie, during 
initial epoch (w = 1), in a region to the right of the second positive root. 

We now study a configuration in which interesting possibilities for bounce and collapse arise 
when b{r) < ; Pe2f{r) > 0. 

If the discriminant of the quintic polynomial is positive, then from the Descartes' rule of signs, 
there are two positive roots 71 (r) and 72 (r) [9], and the space of allowed dynamics is [0,71] and 
[72, 00). The region (71, 72) is forbidden. Shells in the [0, 71] region initially, always become singular. 
Shells initially belonging to the region [72,00) will undergo a bounce and subsequent expansion, 
starting from initial collapse. This bounce occurs when their geometric radius approaches Rbounce = 
r72. If the initial data is chosen such that the discriminant vanishes, then the positive roots are 
equal and there is no forbidden region. 

For 6(r) < ; P6»2/(^) > a shell can evolve in three different possible ways: 

i) If the effective potential for a shell has two positive roots, i. e. if it crosses the x axis twice 
in the range [0, 1] (see FIG.l), then the shell bounces off, as seen by differentiating equation ([13]) 




(14) 




(15) 



to get 2vv + v'^ = —V,y. Since near the turning point V,y < and v > 0, we have v > 0. 



ii) If V{r, u) < in the whole range [0, 1], then the shell will reach the singularity aX v = 0. 

iii) When the potential has double roots in [0, 1], then it indicates that the shell is in critical 
collapse condition. 



In order to understand the possible dynamical evolutions clearly, we now consider one specific 
configuration of the initial data, in which po^ and b{r) are kept fixed and mo is allowed to vary. 
We take 6(r) = — |6o|, then close to the center effective potential can be written as 



V 



P^r\5 ^ _ bor^)v^ - bov + mo(l - m2r'^) 

3mo 



(16) 



We will see that for the above effective potential, depending on the values of initial parameters, 
three different types of evolutions of the collapsing cloud are possible in general. 



Case A 



This type of evolution is depicted in FIG. 2 which shows the effective potential for the central 
shell (r = 0) and two other shells. Initially we have v = 1, v < and as the collapse proceeds the 
value of V decreases, reaching a minimum where potential becomes zero. Then it again increases 
so that V > 0. The effective potential moves towards positive side as one goes away from the 
center and the potential has two positive roots for all the shells. There is complete bounce of 
collapsing shells and Minkowski space is left behind. 



Case B 



As mo is increased further, we reach a value at which the potential for the central shell just 
touches the x axis. We call this value as the critical value moc of the initial parameter mo (for the 
other two parameters fixed). This is shown in FIG. 3. The outer shells still have positive potential 
which indicates forbidden region, therefore, all those shells will bounce back. This configuration 
is the boundary point between the complete dispersal and black hole formation. 



Case C 
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FIG. 2: Effective potential for different comoving radii in Case A type configuration. Here rriQ 
pgrj — 4.5, 6o = —3.0 and m2 = O.I. 
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FIG. 3: Effective potential for different comoving radii for Case B type configuration. Here mo = 1.3333, 
P02 = 4.5, bo = —3.0 and m2 = O.I. 



If we increase rriQ even further, the potential for the central shell (FIG. 4) is negative for the 
whole range of v which will allow the central shell to reach the singularity at v = 0. If we increase 
r, the potential minima goes up, and there is a value of r at which the effective potential just 
touches the x axis. We call this radius the critical radius r^ of the collapsing cloud for the chosen 
set of initial numbers. Now all the shells from r = Q to r = r^ will reach the singularity and will 
contribute to the mass of the black hole formed, but the shells with comoving coordinate more 
than Tc bounce off. We numerically confirm this by integrating equation (jlSp . Physical radius for 
various shells is plotted in FIG. 5. 

We can understand qualitatively why the solution for Case 5 is a massless naked singularity. 



FIG. 4: Effective potential for different comoving radii in Case C type configuration. Here mo = 1.40, 
— 4.5, bo — —3.0 and m2 = 0.1. 
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FIG. 5: Physical radius of different shells in Case C type evolution, where mo = 1.35, |6o| = 3.0, m2 = 0.1 
and pe^ ~ 4.5. Shells with comoving radius less than the critical radius reach singularity, while those with 
larger than critical value escape to infinity. 

Suppose we are in initial data parameter range of Case C regime and we decrease the value of ruQ. 
As rriQ approaches ttiqc fewer and fewer shells reach singularity. Thus if we approach the threshold 
of black hole formation from the initial data of Case C region, we can see that when mg is very 
very close to moc, almost only the central shell reaches singularity and other shells are stopped 
from collapsing and bounce back at the value of v which is the root of the potential for those shells. 
Apparent horizon is given by equation F = R, which implies Vah ~ f'^ for the shells near the 
center. A shell is trapped if it reaches value of v smaller than Vah- Now, as all outer shells reach 
a minimum value vo{r) >> Vah, those are not trapped. Only the central shell is singular when 
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mo is nearly equal to moc and trapped surfaces are not formed. As F = r^A^(r), the mass which 
contributes towards singularity formation is zero which means the singularity is massless and also 
as the trapped surfaces are not formed, the central singularity is visible. 

We note that in the collapsing cloud if a shell bounces then all the shells with larger value of 
comoving radius will also bounce, i.e. there is no shell-crossing. To see this, we can integrate near 
the center numerically to obtain R' as well as all the metric functions. We see that R' remains 
positive during the collapse evolution which is plotted for one particular initial data set in FIG. 6. 
We find all solutions for various values of initial data and it was seen that for the critical solution 
as well as for solutions near the critical point shell cross does not occur. 




r=0.05 
r=0.10 
r=0.15 
r=0.20 



FIG. 6: R' for different shells in Case C type evolution, where mg — 1.35, \bo\ = 3.0, m2 — 0.1 and 



4.5. 



The expression for critical radius in terms of the initial parameters can be obtained from the 
condition that at critical radius, effective potential just touches the e. the quintic poly- 

nomial has double roots. For a quintic polynomial a^x^ + a^x^ + aix + ao = , the discriminant 
is given as follows [lo| . 



A = as [ao(108a| - dOOaiala^ + 2000a?a3a^ + 3125agai) + 16a? - 128afaia5 + 256afa: 



5„21 
5 



(17) 

The condition that A vanishes at the double root gives, from equation (jl6p . the expression for 
critical radius as. 



= mo -mo, 

" mo(0.833&o + 2m2) + ci(0.226[j - 1.04moPe2) ' 
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FIG. 7: Figure represents behavior of critical radius with the density parameter, keeping other parameters 
values unchanged. 

where ci = (3pg.^ + 10m2) and moc = Behavior of the critical radius with the parameter 

mo can be seen in FIG. 7. In the model considered here, the mass function, which is twice the 
Misner-Sharp mass, depends only on r. The total mass which collapsed to form singularity from 
the regular initial profile is ^rj?moe~™^''^. The mass of black hole would then be 



where Cm is a constant near the threshold. We can fix any of the two parameters mo, 6o, and pg^ 
and vary the third one to obtain the expression for critical radius. Following the same procedure 
as earlier, it is easily seen that for bo and pg^ the same scaling relation and exponent exists near 
the threshold of the black hole formation. Therefore, in general, for a parameter ^, we can write 



We like to point out the similarities in the behaviour of this model with the critical behavior in 

n f] 

scalar field gravitational collapse [11], and other matter models [5[: 

i) The three types of evolutions discussed in Case ^-Care very much analogous to the subcritical, 
Critical and super- critical evolutions studied in the above models. In case A the whole cloud 
disperse off and in Case C black hole forms whose mass can be controlled by an initial parameter. 
Case B is the boundary between the case A and Case C configuration which is like the critical 
solution in the above cases. 




(19) 



Mbh = C(.\C - Cc\^ . 



V. DISCUSSION AND CONCLUSIONS 
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ii) A power law behaviour is seen for the mass of the black hole near the threshold of its 
formation. Infinitesimal mass black holes can be formed in this model if the initial parameters are 
finely tuned. 

iii) The same value of exponent is observed for different initial parameters which is like 
universality of the critical solution. 



Our purpose here is to investigate how the initial parameters determine the evolution and the 
end state of collapse. The metric function i^, in spite being a restrictive choice, makes the model 
tractable. The other class of tangential pressure models, i.e non-static Einstein cluster has a 
different form of the function v in comoving coordinates and it can not be written as function 
of R only. We have observed similar behavior at threshold of black hole formation in Einstein 



cluster as well and the same value of the exponent was seen _12|]. The effective potential method 
'or calculating the critical exponent applies to this case as well as to non-static Einstein cluster 



12l |. This suggests that the method applies at least for all mass-conserving systems. However, 
the system considered here has a limitation, and is different from the models considered so far 
for critical behaviour in one respect that it has no radial stress which enables shells to interact 
directly with each other. It will be interesting to explore whether the method we used here could 
be applied for studying such behaviour in the models with a non-vanishing radial pressure as well. 
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APPENDIX A : BLACK HOLE FORMATION FOR COLLAPSE IN CASE C 

CONFIGURATION 



We can rewrite equation (|T2|) as 



= -^e^^^vhQ e^'^o (v^hirv) + M), (20) 

where 

^(^) = R2 ■ (21) 

Integrating the above equation, we get 

tiv, r) = C (22) 
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The time of singularity for a shell at a comoving coordinate radius r is the time when the physical 
radius R{r, t) becomes zero. The shells collapse consecutively, that is one after the other to the 
center as there are no shell-crossings. Taylor expanding the above function around r = 0, we get 



t{v,r) = t{v,0) + r 



dt{v, r) 



dr 

Let us denote 

Xn{v} = 



+ 777 



r^ dH(v,r) 



(23) 

r=0 



(24) 

r=0 



dr"' 

As we have taken the initial data with only even powers of r, the first derivatives of the functions 
appearing in above equations vanish at r = 0, hence we have 

Xi{v) = (25) 



Now we can express the next coefHcient <%2 as 

9ml 0^'"'^ + ^oov + M)3 



^/vdvlPv^ + 2moPe2V^ + 9m^(6o2^^ + "12)] 
^2(v) = - ,. , ^ TZ-i (^oj 

Jv 



where 

/? = 6pg^ + 24mo6o2P6i2 - 3mo/92P6i2 (27) 

We need to determine now whether it is possible to have families of future directed outgoing null 
geodesies coming out of the singularity. In the case when such families do exist which terminate 
in the past at the singularity, and which could reach outside observers, then the singularity will 
be visible. In the case otherwise it is hidden within the black hole. Another way to look at this is 
through the apparent horizon and formation of trapped surfaces in the spacetime. As the collapse 
evolves, if the trapped surfaces form well in advance to the formation of the singularity, then the 
same will be covered. On the other hand, if the trapped surface formation is sufficiently delayed 
during the collapse then the singularity may be naked. The apparent horizon within the collapsing 
cloud is given by the equation , R/F = 1, which gives the boundary of the trapped surface region 
of the space-time. If the neighborhood of the center gets trapped earlier than the singularity, then 
it is covered, otherwise it is naked with non-spacelike future directed trajectories escaping from it. 

In order to consider the possibility of existence of such families, and to examine the nature of 
the singularity occurring at i? = 0, r = in this model, let us consider the outgoing null geodesic 
equation which is given by 

^ = e^-" (28) 
dr 
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We now use here a method which is similar to that given in [13]. The singularity curve is given 
by v{ts{r),r) = 0, which corresponds to R{ts{r),r) = 0. Therefore, if we have any future directed 
outgoing null geodesies terminating in the past at the singularity, we must have ii ^ as t — > 
along the same. Now writing equation (j28p in terms of variables {u = r^jR), we have 



du a 



(29) 



Now in order to get tangent to the null geodesic in the {R, u) plane, we choose a particular 
value of Q such that the geodesic equation is expressed only in terms of (^) . A specific value of 
alpha is to be chosen which enables us to calculate the proper limits at the central singularity. For 
example, for A'i(O) 7^ case, we can choose a = | and using equation ([5]), (and considering that 
i? < 0) , we get 

In the tangential pressure collapse model discussed in the previous section we have 'Vi(O) = 0, and 
hence we choose a = | so that when in limit r ^ 0, t ^ we get the value of tangent to null 
geodesic in the {R, u) plane 




- ^ ^ (31) 

G{^/G+^/H) 

Now note that for any point with r > on the singularity curve ts{r), we have i? — > whereas 
F (interpreted as mass of the object within the comoving radius r) tends to a finite positive value 
once the energy conditions are satisfied. Under the situation, the term F/R diverges in the above 
equation, and all such points on the singularity curve will be covered as there will be no outgoing 
null geodesies from such points. 

Hence we need to examine the central singularity at r = 0, i? = to determine if it is visible or 
not. That is, we need to determine if there are any solutions existing to the outgoing null geodesies 
equation, which terminate in the past at the singularity and in future go to a faraway observer, 
and if so under what conditions these exist. Note that if any outgoing null geodesies terminate at 
the singularity in the past then along the same, in the limit as r ^ 0, t — > we then have from 
equation (26) i? = 0, therefore H = and G = 1 in this limit as F/R vanishes. Let now xq be the 
tangent to the null geodesies in (i?, u) plane, at the central singularity, then it is given by, 

R dR 



Xq = lim lim 

t^ts r^O U du 



(32) 
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Using equation ([3T]) . we get, 

4 = IVm'oMO) (33) 

In the {R, u) plane, the nuh geodesic equation will be, 

R = xqu (34) 

It follows that if ^^2(0) > 0, then that implies that xq > 0, and we then have radially outgoing 
null geodesies coming out from the singularity, making the central singularity to be a visible one. 
On the other hand, if ^^2(0) < 0, we will have a black hole solution. Now for the initial data in the 
range of Case C, we calculate ^^2(0) and find that it remains negative which implies that a black 
hole is necessarily formed. 
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